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ABSTRACT 

We use thermochemical equilibrium calculations to model iron, magnesium, 
and silicon chemistry in the atmospheres of giant planets, brown dwarfs, ex- 
trasolar giant planets (EGPs), and low- mass stars. The behavior of individual 
Fe-, Mg-, and Si-bearing gases and condensates is determined as a function of 
temperature, pressure, and metallicity. Our results are thus independent of any 
particular model atmosphere. The condensation of Fe metal strongly affects 
iron chemistry by efficiently removing Fe-bearing species from the gas phase. 
Monatomic Fe is the most abundant Fe-bearing gas throughout the atmospheres 
of EGPs and L dwarfs and in the deep atmospheres of giant planets and T 
dwarfs. Mg- and Si-bearing gases are effectively removed from the atmosphere 
by forsterite (Mg2Si04) and enstatite (MgSiOs) cloud formation. Monatomic 
Mg is the dominant magnesium gas throughout the atmospheres of EGPs and 
L dwarfs and in the deep atmospheres of giant planets and T dwarfs. Silicon 
monoxide (SiO) is the most abundant Si-bearing gas in the deep atmospheres 
of brown dwarfs and EGPs, whereas SiH4 is dominant in the deep atmosphere 
of Jupiter and other gas giant planets. Several other Fe-, Mg-, and Si-bearing 
gases become increasingly important with decreasing effective temperature. In 
principle, a number of Fe, Mg, and Si gases are potential tracers of weather or 
diagnostic of temperature in substellar atmospheres. 



Subject headings: astrochemistry — planets and satellites: individual (Jupiter) 
— stars: low-mass, brown dwarfs — stars: individual (Gliese 229B, HD 209458) 
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Introduction 



Lodders fc Fegley 



Ther mo chemical mo dels have been used to describe alkali (|Lodderslll999l). ti tanium and 
vanadium (lLoddersll2002[). carbon , nitrogen, and oxv een (ILodders fc Feglevll2002l) . sulfur and 



2006 



Loddersll2009r ) 



phosphorus ( iVisscher et al.ll2006l ). condensate (e.g., 
chemistry in the atmospheres of gas giant planets, brown dwarfs, and low-mass dwarf stars. 
Here we continue and extend these previous studies by using thermochemical equilibrium cal- 
culations to model the chemical behavior of Fe, Mg, and Si in substellar objects. Iron, mag- 
nesium, and silicon are the most abundant rock-forming elements in a solar composition gas, 
and - condensed as iron metal (Fe), and forsterite (Mg2Si04) and enstatite (MgSiOs) - will 
produce the most massive cloud layers in substellar atmospheres. Cloud formation strongly 
affects the optical and infrared spectra of substellar objects by removing gases from the over- 



lying atmosph e re and by intro duciiig solid or liquid cloud pa r ticles fe.g.. iMar 



LoddersI Il999[ iBurrows et al.l I2OOOI : iLodders fc Fegleyl I2OO6I : IVisscher et al. 
2009h . 



ev et al 



2006 



1996 



Lodders 



The large and growing number of discovered brown dwarfs (>750) and extrasolar plan- 
ets (>400) makes it imp ractical to model the thermochemist ry of all objects individ ually 
( iLodders fc FeglevI 120021). as was don e for Jupiter and Saturn (IFegley fc LoddersI 119941 ) and 



Gliese 229B flFeg. 



ev fc Lodderslll996l). Instead, we adopt an a pproach similar to that of our 



previous papers ( iLodders fc Fegleyll2002t IVisscher et al.ll2006[ ) and determine the abundance 



of each chemical species as a function of pressure, temperature, and metallicity, plotted in 
abundance contour diagrams. Our results are thus independent of any particular pressure- 
temperature profile, and in principle, the atmospheric profile for any object may be super- 
imposed on the abundance diagrams to determine its equilibrium atmospheric chemistry. In 
some instances, the behavior of key gases may be diagnostic of atmospheric temperature 
and/or tracers of weather in substellar objects. 

We begin with a brief description of our computational method ([|2]), and then present 
our results for iron chemistry in substellar atmospheres in fJ21 We first give an overview of iron 
chemistry in a solar composition gas and identify important gases and condensates ( §3.ip . 
This is followed by more detailed discussion of the chemical behavior individual Fe-bearing 
gases as a function of temperature, pressure, and metallicity ( §3.2p . Wherever possible, we 
note relevant spectroscopic observations of Fe-bearing gases in substellar objects. Our results 
for iron are then summarized by illustrating Fe gas chemistry along the atmospheric profiles 
of representative substellar objects ( §3.3p . A similar approach to magnesium and silicon 
chemistry follows in f|H and ^ respectively. We conclude with a brief summary in ^ 
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2. Computational Method 



Thermochemical equilibrium calculations were performed using a Gibbs free energy min- 
imiza tion code, previously used for modeling the atmospheric chemistr y of Saturn (IVisscher et al 



20061 ) and sulfur and phosphorus chemistry in substellar objects (IVisscher et al.l |2006| ). 



Thermodynaraic data for t he equilibrium calculations were taken from the compilations of 
Gurvich et al.Nl989-1994f l. lRobie fc Hemingway! Jl995h . the fourth edition of the JANAF Ta- 
bles ( 1Chaselll998r). and the thermodyn amic database maintained in the Planetary Chemistry 
Laboratory (ILodders fc Fegleyl |2002| ). This database includes additional thermodynamic 
data from the literature for compounds absent from the other compilations. 



All calculations were conducted using elemental abundances from iLodderd (120031 ) for a 
solar system (i.e., protosolar) composition gas. The effect of metallicity on sulfur and phos- 
phorus chemistry was examined by running computations at [Fe/H] = -0.5 dex (subsolar), 
[Fe/H] = dex (solar), and [Fe/H] = +0.5 dex (enhanced) metallicities. The metallicity fac- 
tor, m, is defined as logm = [Fe/H]. We assume that the elemental abundance ratios for Mg, 
Si, and other elements of interest vary similarly wit h [Fe/H] (e.g., [Mg/H] [Si/H] !^ [Fe/H]) 



over the range of m etallicities considered here (see lEdvardsson et al. 



1993; Chen et al. 2000 



Huang et al.l l2005l ) . When considering the chemical behavior of individual gases, we focus 
on higher temperatures (800 K and higher), where thermochemical processes are expected to 
dor ninate over disequilib rium processes such as photochemistry or atmospheric mixing (e.g.. 



see 



Visscher et al.ll2006f ) 



We assume that condensates settle gravitationally into a cloud layer and are removed 
from the cooler, overlying atmosphere. This equilibrium cloud condensate scenario for the 
deep atmosphere s of giant planets and brown dwarfs is supported by several lines of ev- 



idence (e.g., see iLodders fc Fegleyl l2006l . and references therein). First, the presence of 



germane (GeH4) and the absence of silane (SiH4) in the upper atmospheres of Jupiter and 
Saturn (even though Si is expected to be much more abundant than Ge) can be explained 
by the removal of Si from th e gas into silicate clouds deeper in the atmo sphere, whereas 
Ge remains in the gas phase (IFegley fc PrinnI Il988l : iFegley fc Lodderd [l99J) . Secondly, the 
detection of H2S in Jupiter's troposphere by the Galileo entry probe indicates that Fe must 
be sequestered into a cloud layer at deep atmospheric levels, beca use the formation of FeS 
would otherwise rem ove H2S from the gas above the ~ 700 K level ( ILodders &: Fegleyl 12002 



Visscher et al 



20061) . Third, absorption from monatomic K gas in the spectra of T dwarfs 



( jBurrows et al. 



2OOOI : iGeballe et al.ll200ll ) requires the removal of Al and Si at deeper atmo- 



spheric levels, because K would otherwi se be removed frorn the observable atmosphere by the 
condensation of orthoclas e (KAlSiciO«) (ILodders fc Feglevll2006r). The p r esence of monatomic 



Na gas in brown dwarfs (IKirkpatrick et al 



1999 



Burgasser et al.l l2003l : McLean et al.l 12003 
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Cushing et al.ll2005l ) also suggests Al and Si removal, because albite (NaAlSisOs) condensa- 



tion would otherwise effectively remove Na from the observable atmosphere. Furthermore, 
the removal of Na by Na2S cloud formation is consistent with the observed weakening of 
Na atomic lines throughout the L dwarf spectral sequence a r id their disappearan ce in early 
T dwarfs (e.g., iKirkpatrick et al.lll999[ iMcLean et aD boosl : Ivisscher et aD booeh . Finally, 
as we note below, the disappearance of iron, magnesium, and silicon spectral features in 
later spectral types is consistent with removal of these elements into cloud layers. In our 
thermo chemical model, the abundances of Fe-, Mg-, and Si-bearing gases above the clouds 
are computed assuming saturation (equilibrium) vapor pressure. If supersaturation occurs, 
a condensate will form and settle toward the cloud layer to restore equilibrium. 



3. Iron Chemistry 



3.1. Overview of Iron Chemistry 

Figure [T]illu strates model atmospheric profiles for a n M dwarf (T^ff = 260 K, log g = 5.0; 



Tsuji et al. 



1996h . an L dwarf (Teg = 1800 K, l ogo = 5.0:[Burrows et al.l booeh. the hot, close 



in (Pegasi ) planet HD209458b {T^s = 1350 K; liro et al.ll2005f ). the T dwarf Gliese 229B {T^s 
= 960 K; iMarley et al.lll996l ). and Jupiter {Tes = 124 K), indicated by dashed lines. We 
note that Jovian atmospheric chemistry differs slightly than that for a solar-metall icity gas 
because Jupiter has a heavy e lement enrichment comparable to [Fe/H] ^ +0.5 dex (ILodders 
I999I : iLodders fc Feglevl[2002l ). 



Also shown in Figure [T] are lines indicating where A(CH4) = A(CO) and A(H2) = A(II). 
These boundaries are important because carbon and oxygen affect the chemical behavior 
of many Fe, Mg, and Si-bearing gases. Methane is the dominant carbon-bearing gas in 
Jupiter and T dwarfs (such as Gliese 229B) whereas CO is the dominant carbon-bearing 
gas in L dwarfs and Pegasi planets (such as HD209458b). Molecular hydrogen dissociates 
into monatomic H at high temperatures and low pressures (lower right corner of Figured]). 
However, H2 is the dominant form of hydrogen in substellar objects, and we therefore take 

~ -^SH ~ 0.84 throughout the following. This approximation holds for metallicities up 
to [Fe/H] +0.5 dex; at higher metallicities the H2 mole fraction abundance decreases as 
the relative abundance of heavy elements increases (e.g., ~ 0.82 at [Fe/H] ^ +1.0 dex). 

The dotted line in Figure [1] shows the condensation curve for Fe metal, with an open 
circle denoting its normal melting point (1809 K). Iron condensation occurs via 



Fe = Fe(s,l). 



(1) 
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The equilibrium condensation temperature (Tcond) of Fe is approximated by 

lOVTcond(Fe) = 5.44 - 0.48 log Pt - 0.48[Fe/H], (2) 

where higher pressures and/or metallicities lead to higher Fe condensation temperatures. 
Iron cloud formation effectively removes nearly all iron from the atmosphere, and the abun- 
dances of Fe-bearing gases above the clouds rapidly decrease with altitude. 

Monatomic Fe gas is the dominant Fe-bearing gas in a protosolar composition gas until 
it is replaced by Fe(0H)2 at low temperatures and high pressures. The conversion between 
Fe and Fe(0H)2 is represented by the net thermochemical reaction 

Fe + 2H2O = Fe(0H)2 + H2, (3) 

and the solid line dividing the Fe and Fe(0H)2 fields in Figure [T] indicates where these gases 
have equal abundances A(Fe) = A(Fe(0H)2) ~ |EFegas, where SFcgas is the total amount 
of iron in the gas. The position of this line is given by 

log Pt = 11-94 - 12088/r-2[Fe/H], (4) 

showing that an increase in metallicity will shift the Fe-Fe(0H)2 boundary to higher temper- 
atures and lower pressures. In other words, an increase in metallicity increases the stability 
field of Fe(0H)2 for otherwise constant conditions. 



3.2. Chemical Behavior of Iron-bearing Gases 

3.2.1. Monatomic Iron, Fe 

The mole fraction abundance of monatomic Fe as a function of pressure and temperature 
is shown in Figure [2^. As mentioned above, monatomic Fe is the dominant Fe-bearing gas 
over a wide range of pressures and temperatures in a solar system composition gas. Below 
the Fe clouds, the abundance of Fe gas is representative of the total iron content of the 
atmosphere {X-pc ~ -^efg), and its abundance is given by 

logXFe~-4.24+[Fe/H]. (5) 

Upon Fe metal condensation, the amount of iron in the gas rapidly decreases and the equilib- 
rium Fe gas abundance is governed by its vapor pressure over solid or liquid iron, represented 
in reaction ([1]). The mole fraction abundance of monatomic Fe above the clouds is given by 



logXpc ~ 7.23 - 20995/T - log Pt, 



(6) 
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inversely proportional to Pt- The Fe gas abundance is independent of metallicity in this 
region because it depends solely on the temperature-dependent vapor pressure of iron. As 
described below (see §3.2.21 and §3.2.3p . this expression can be used with chemical equilibria 
to determine the abundances of other iron gases in substellar atmospheres. Equations giving 
the abundance of Fe (and other Fe-bearing gases) as a function of temperature, pressure, 
and metallicity below and above the Fe clouds are listed in Table [H 

Neutral monatomic Fe possesses several hundred spectral lines in the J and H bands 
(IMelendez fc Barbuyl 119991 ) , an d Fe I features are ob served in the spectra of brown dwarfs 



and low-mass dwarf stars (e.g. ICushing et al.ll2005l ). The strong Fe feature at 1.189 /im 



weakens in mid- to late-type M dwarf spe ctra and generally disappears in mid-type L dwarfs 
(IMcLean et al.ll2003l : ICushing et al.ll2005l ). This trend is consistent with the removal of iron 
from the gas into an Fe metal cloud deck located deeper be low the photosphere in objects 
with low effective temperatures (e.g. iBurgasser et al.ll2002bl ). 



3.2.2. Iron Hydride, FeH 

The chemical behavior of FeH in a protosolar composition gas is illustrated in Figure 
The conversion between Fe and FeH occurs via the net thermochemical reaction 

Fe + O.5H2 = FeH (7) 

Here we show how chemical equilibria may used to derive equations giving the abundance of 
Fe-bearing species as a function of temperature, pressure, and metallicity. Rearranging the 
equilibrium constant expression for reaction ([7]) yields 

log XpeH = log Xpe + 0.5 log Xh, +0.5 log Pt + log (8) 

Using the Fe abundance from equation ([5]), the temperature dependence of (logiij^ ~ 
-1.85 - 1905/T from 800 to 2500 K), and the hydrogen abundance (Xr^ ^ 0.84), the FeH 
abundance between the H2-H boundary and the Fe cloud deck is given by 

logXpeH ^ -6.13 - 1905/r + 0.51ogPr + [Fe/H], (9) 

proportional to P^'^ and m. At high temperatures (> 2000 K) and low pressures (< 10^^ bar) 
as monatomic H becomes increasingly abundant the H2 abundance begins to decrease near 
the H2-H boundary. The reduced H2 abundance, in turn, reduces FeH formation via reaction 
([7]) and changes the chemical behavior of FeH at high temperatures and low pressures (i.e., 
lower right corner of Figure [2^) . 
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Above the iron clouds, the Fe abundance in equation ([6]) is used in equation ([8]) to give 



logXpcH ~ 5.34 - 22900/T - 0.5 log Pt, 



(10) 



where XpeH is proportional to ' and independent of metallicity. The differences in 
chemical behavior of FeH below (equation |9]) and above (equation [10]) the Fe cloud deck are 
illustrated in the shape of the FeH abundance contours in Figure Wp- The sharp bends in 
the contours correspond to the Fe condensation curve in Figure [H 

Absorption bands from FeH are common in the near-infrared spectra of brown dwarfs 



(e.g. iKirkpatrick et al.l Il999l : McLean et all 120081 : ICushing et al.l 120031 |2005|). The most 
prominent FeH feature is the band located at .9896 /im, which weakens throughout the 
L dwarf spectral seq uence and in early T dwarfs (IMcLean et al.ll2003l : iNakajima et al.ll2004j : 
Gushing et al.ll2005l ). consistent with the removal of iron into a cloud located deeper and 
deeper i n the atmosphere. This fe ature unexpectedly strengthens again in ra id-type T dwarf 
spectra (IBurgasser et al.ll2002al Jbl: INakajima et al.ll2004j : ICushing et al.ll2005l ). prompting dif- 
ferent explanations for the shape of the color- magnitude diagram for brown dwarfs near 
the L-T transition. INakajima et al.l ( 12004 ) suggested that the strengthening FeH bands are 
caused by upward convective mixing of FeH ga s from d eeper levels where it is more abundant. 
Howev er, as pointed out by IBurgasser et al.l ( l2002bl ) and demonstrated iLodders fc Fegley 
( 120061 ). the fragile Fe — H bond is unlikely to survive convective upwelling in a T dwarf atmo- 
sphere. Instead, the observations are plausibly explained by cloud disruption and clearing 
(in a ~1 /xm window) which allows the observa tion of FeH gas at deep atmospheric levels 
dSurgasser et al.ll2002b[ Ihodders fc Fedevl 120061 ). 



3.2.3. Iron Monohydroxide, FeOH 

Mole fraction contours of FeOH are illustrated in Figure EJ:;. The equilibrium between 
Fe and FeOH is represented by the net thermochemical reaction 

Fe + H2O = FeOH + O.5H2, (11) 

and expressions giving the FeOH abundance as a function of temperature, pressure, and 
metallicity are listed in Table [1] Above the iron clouds, the curvature in the FeOH abundance 
contours along the CH4-CO boun dary results from the effe ct of carbon chemistry on the H2O 



abundance in reaction ( ITTl) (e.g. iLodders fc Fegley! l2002l ). For example, the atmospheric 



water abundance may be written as 

logXH,o = logX^j^o + [Fe/H], (12) 



- 8 - 



where X{^^q is the water abundance in a solar-metalhcity gas, which is logj^gX^^o ~ 10~^'^^ 
inside the CH4 field and X^^^q ^ ]^q-3.58 j^side the CO field in Figure [H At temperatures 
and pressures near the CH4-CO boundary, X^^q may be derived from CH4-CO equilibria 



(iLodders fc Fegleyll2002l ). Rearranging the equilibrium constant expression for reaction flTTl) 
gives 

logXpeOH = logXpe + logXn^o " O.SlogXna + 0.5 log Pt + log%i} (13) 
Substituting for the Fe abundance from equation ([HD, the H2O abundance from equation 
(fT2|) . the H2 abundance {Xu^ ^ 0.8384), and the temperature dependence of I^YTl (^^S ^^TTl ~ 
-2.25 + 2468/T from 800 to 2500 K), equation ([13]) becomes 

logXpeOH ~ 5.02 - 18527/T - 0.5 log + 2[Fe/H] + logX^^Q. (14) 

This expression gives the FeOH abundance above the iron clouds and includes the effect of 
metallicity on the atmospheric water abundance. The kinks in the FeOH contours denote 
the position of the Fe condensation curve. As demonstrated in equation f[T^ and shown in 
Figure [2]3, the mole fraction abundance of FeOH is proportional to P^°'^ throughout this 
region. 



3.2.4. Iron Hydroxide, Fe(0H)2 

The chemical behavior of Fe(0H)2 as a function of temperature and pressure is illus- 
trated in Figure [2]i. The equilibrium conversion between Fe and Fe(0H)2 is represented by 
reaction ([3]). Below the Fe clouds (kinks in abundance contours), the Fe(0H)2 abundance 
increases with total pressure, and has a strong dependence on metallicity (see Tabled]). 

The Fe(0H)2 abundance in reaction is very sensitive to the H2O abundance since 
the formation of one mole of iron hydroxide requires two moles of water. As a result, there is 
a large shift in the Fe(0H)2 abundance contours above the clouds (inflections in Figure [2]1) 
when moving between the CH4 and CO flelds. Within each fleld, the Fe(0H)2 abundance 
above the clouds is essentially pressure-independent (see Tabl^ and is therefore, in principle, 
diagnostic of atmospheric temperature. 

As shown in Figure [H Fe(0H)2 becomes the most abundant Fe-bearing gas at low 
temperatures and high pressures in a solar composition gas. Iron hydroxide is therefore 
expected to be the most abundant iron gas at T < 1600 on Jupiter ([Fe/H] f» -1-0.5) and 
T < 1070 on Gliese 229B. However, we emphasize that the amount of iron remaining in the 
gas is greatly diminished by the condensation of Fe metal at deeper atmospheric levels. For 
example, even where it is the most abundant Fe-bearing gas, the predicted Fe(0H)2 mole 
fraction is ~ 10~^^'^ at the 1000 K level in Jupiter's atmosphere. 
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3.2.5. Iron Monoxide, FeO 

Mole fraction abundance contours for FeO are shown in Figure [3^ for a solar composition 
gas. Monatomic iron reacts with water to form FeO via the net reaction 

Fe + H2O = FeO + H2. (15) 

Between the H2-H boundary and the Fe cloud deck, Xpeo is proportional to and is inde- 
pendent of pressure (see Tabled]). In principle, the FeO abundance in this region is therefore 
diagnostic of either temperature or metallicity. Above the clouds, the FeO abundance de- 
creases with pressure and increases with metallicity. As for FeOH and Fe(0H)2, the FeO 
abundance is sensitive to the H2O abundance and the FeO abundance contours display a 
subtle shift along the CH4-CO equal abundance boundary. 



3.2.6. Iron Sulfide, FeS 



The abundance of FeS g function of temperature and pressure is illustrated in 

Figure The chemical behavior of FeS is governed by the reaction between Fe and H2S: 



Fe + H2S = FeS + H2. 



(16) 



Above the iron clouds, H2S is the dominant sulfur-bearing gas (logXnaS ~ —4.52 + [Fe/H]), 
and the FeS abundance is proportional to P^^ and m (see Tabled}. Below the Fe cloud deck, 
the FeS abundance is ~ 10 ppb and decreases at high er temperatures and lower pressures 
as H2S is replaced by SH and S ( IVisscher et al.ll2006l ). Iron sulfide is predicted to be the 
second or third most abundant iron gas throughout L dwarf atmospheres and in the upper 
atmospheres (< 1500 K) of Pegasi planets. 



3.3. Iron Chemistry in Substellar Objects 

Figures H^-HJi summarize the iron gas chemistry in four representative substellar objects: 
Jupiter, the T dwarf Gliese 229B, the Pegasi planet HD209458b, and an L dwarf (Teg = 1800 
K). Iron chemistry in substellar objects is strongly affected by Fe metal condensation at deep 
atmospheric levels. Monatomic Fe is the dominant Fe-bearing gas in the deep atmospheres of 
Jupiter and Gliese 229B, and throughout the atmospheres of HD209458b and the L dwarf. A 
number of Fe-bearing gases become relatively more abundant in objects with lower effective 
temperatures. On Jupiter and Gliese 229B, Fe(0H)2 replaces Fe at lower temperatures. The 
Fe(0H)2 and FeO abundances are pressure-independent and thus potentially diagnostic of 
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atmospheric temperature, respectively, above and below the iron clouds. For all four objects, 
iron hydride (FeH) is the second most abundant Fe-bearing gas at deep atmospheric levels, 
until it is replaced at lower temperatures by Fe(0H)2, FeOH, or FeS. Because of its strong 
absorption and relatively high a bundance near the iron condensation level, FeH is a tracer 
of weather in brown dwarfs (e.g. iBurgasser et al.l l2002bl ) . Other Fe-bearing gases (e.g. Fe, 
FeOH, FeS) are potential tracers of weather since they typically achieve their maximum 
abundance near iron cloud base. 



4. Magnesium Chemistry 

4.1. Overview of Magnesium Chemistry 

Figure gives an overview of Mg chemistry as a function of pressure and temperature 
in a protosolar composition gas. The dotted lines labeled Mg2Si04(s,l), MgSi03(s,l), and 
MgO(s) show the equilibrium condensation curves for forsterite, enstatite, and periclase, 
respectively, and the open circles denote the normal melting temperatures for forsterite 
(2163 K) and enstatite (1830 K). Forsterite (Mg2Si04) condenses via the net thermochemical 
reaction 

2Mg + 3H2O + SiO = Mg2Si04(s,l) + 3H2, (17) 

and its condensation temperature as a function of pressure and metallicity is approximated 
by 

10V^cond(Mg2SiO4) ^ 5.89 - 0.37 log Pt - 0.73[Fe/H], (18) 
At slightly lower temperatures, enstatite (MgSiOa) condensation occurs via the net reaction 

Mg + 2H2O + SiO = MgSi03(s,l) + 2H2. (19) 

The enstatite condensation curve is approximated by 

10V^cond(MgSiO3) ^ 6.26 - 0.35 log Pt - 0.70[Fe/H], (20) 

where higher pressures and/or metallicities lead to higher condensation temperatures. At 
very high pressures, forsterite and enstatite condensation temperature are depressed as SiO is 
replaced by SiH4 (see §5.ip . Periclase (MgO) condenses via the net thermochemical reaction 

Mg + H20 = MgO(s) + H2, (21) 

at pressures greater than 10^'^ bar. The condensation curve for akermanite (Ca2MgSi207) 
is not shown here because most Ca is expected to be removed at deeper atmospheric levels 
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by the condensation of refractory calcium aluminates (e.g., iLoddera |2002| ). If no calcium 
is removed, Ca2MgSi207 condensation would consume ~ 3% of the total atmospheric Mg 
inventory. In the same way, olivine ((Mg,Fe)2Si04) and fayalite (Fe2Si04) are not expected 
in substellar atmospheres because nearly a ll Fe is removed frora the gas phase by iron rn etal 
condensation at higher temperatures (e.g.. ICushing et al.ll2006l : iLodders fc Fegleyl 120061 ) 



Magnesium-silicate cloud formation is very effective at removing nearly all (>99%) 
magnesium from the atmosphere, and the abundances of Mg-bearing gases rapidly decrease 
with altitude above the clouds. This behavior is generally consistent with the disappearance 
of Mg spectral features by early- type L dwarfs (see §4.2.ip . Furthermore, ICushing et al. 
(120061 ) find a Si — O absorption feature at 10 fim in mid-type L dwarfs which is consistent 
with the presence of silicate grains and a weak 9.17 fim feature tentatively attributed to 
crystalline enstatite. 

Monatomic Mg is the dominant Mg-bearing gas in substellar atmospheres until it is 
replaced by Mg(0H)2 at low temperatures and high pressures. The conversion between Mg 
and Mg(0H)2 is represented by the net thermochemical reaction 

Mg + 2H2O = Mg(0H)2 + H2, (22) 

and the solid line in Figure \5\ indicates where Mg and Mg(0H)2 have equal abundances 
A(Mg) = A(Mg(0H)2) ~ |SMggas, where EMggas is the total amount of magnesium in the 
gas. The position of the equal-abundance line is given by 



log Pt = 10.96 - 10267/r - 2[Fe/H], 



(23) 



where an increase in metallicity shifts the Mg-Mg(0H)2 boundary to higher temperatures 
and lower pressures. 



4.2. Chemical Behavior of Magnesium-bearing Gases 

4.2.1. Monatomic Magnesium, Mg 

Figure ^ shows the chemical behavior of monatomic Mg as a function of pressure and 
temperature in a solar-metallicity gas. The chemical behavior of Mg is strongly affected by 
silicate cloud formation. Below the forsterite clouds, the abundance of Mg gas is given by 

logXMg~-4.15 + [Fe/H], (24) 

and comprises nearly 100% of the total elemental Mg content in the atmosphere [Xug ~ 
-^EMg)- Upon forsterite and enstatite condensation, the magnesium gas abundance is gov- 
erned by its vapor pressure over rock, represented by reactions f|T7|) and ([T9l) . Curvature 
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in the Mg mole fraction contours in Figure |6^ occurs along the CH4-CO boundary, which 
affects the H2O abundance in reactions f|T71) and f|T9|) . The Mg abundance is therefore 



logXMg ^ 11.37 - 27250/T - logP^ - [Fe/H]. 
in CH4-dominated objects, and 

logXMg ~ 11.83 - 27250/T - log Pt - [Fe/H]. 



(25) 



(26) 



in CO-dominated objects. We can use these expressions along with chemical equilibria to 
determine the abundances of other Mg-bearing gases in substellar atmospheres. Equations 
giving the abundance of Mg-bearing g function of pressure, temperature, and metal- 

licity in a protosolar composition g listed in Table [2l 

Several Mg absorption lines are pr esent in the near infrared spectra of M dwarfs, includ- 
ing th e prominent feature at 1.183 fim (jjones et al.lll996l : iMcLean et al.ll2003l : ICushing et al. 
2OO5I ). These features wea ken in mid- to late- type M dwarfs and generally disappear by 
~ LI (ICushing et al.l 120051 ). consistent with the removal of magnesium into Mg2Si04 and 
Mg SiO.s clouds lo c ated at increasingly grea ter depths below the observable atmosphere (e.g.. 



see 



Loddersll2004l : iLodders fc Feglevll2006f ) 



4-2.2. Magnesium Hydride, MgH 

Mole fraction contours for MgH are illustrated in Figure (6)3. The MgH abundance is 
governed by equilibrium with monatomic Mg, via the net thermochemical reaction 



Mg + O.5H2 = MgH 



(27) 



Between the H2-H boundary and the Mg- silicate cloud base, the MgH abundance is propor- 
tional to P^'^ and (see Table [2]). In contrast, the MgH abundance above the magnesium- 
silicate clouds is proportional to Pj7°'^ and m~^. The sharp bends in the MgH abundance 
contours occur at the condensation temperature of forsterite. Slight inflections in the MgH 
abundance occur along the CH4-CO boundary, because the Mg abundance in reaction (!27|) 
is sensitive to the H2O abundance in reactio ns (fTTll and (flQll. Magnesium hydride is an im- 



portant opacity source from 0.44 to 0.56 fim (jWeck et al. 



20031), and MgH bands at 0.48 and 



0.52 nm have been found in the optical spectra of bright L dwarfs and extreme sub dwarfs 



flKirkpatrick et al.lll999l : iReid et al.ll2000 



Kirkpatrickll2005f ). 
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4-2.3. Magnesium Monohydroxide, MgOH 

The chemical behavior of MgOH as a function of pressure and temperature is shown in 
Figure [6t. The equihbrium abundance of MgOH is governed by the reaction 

Mg + HaO = MgOH + 0.5H2, (28) 

decreasing with temperature and proportional to P^'^ and below the Mg-silicate clouds. 
The kinks in the MgOH contours occur where forsterite condenses. Above the clouds, XngOH 
is proportional to P^^"'^ (see Table [2]). Interestingly, the MgOH abundance in this region 
is independent of metallicity and shows no shift at the CH4-CO boundary because these 
effects cancel out in reaction (128!) . For examp le, the H^O abundance slightly decreases when 



moving from the CH4 to the CO field (e.g. iLodders fc Fegleyl 120021 ). However, as can be 



seen by comparing equations (12511 and ( ]26l) . there is a corresponding increase in the Mg 
abundance, so the resulting MgOH abundance in reaction (I28p remains unaffected by carbon 
chemistry. Magnesium monohydroxide is typically the third most abundant Mg-bearing gas 
in substellar atmospheres. 



4.2.4. Magnesium Hydroxide, Mg(0H)2 

Mole fraction abundance contours for Mg(0H)2 are shown in Figure Oi for a solar- 
metallicity gas. The equilibrium conversion between Mg and Mg(0H)2 is represented by 
reaction fl22|) . Below the Mg-silicate clouds, the Mg(0H)2 abundance decreases with tem- 
perature and pressure, and has a very strong (m^) dependence on metallicity. 

Above the magnesium-silicate clouds (the kinks in Figure [Hli), the Mg(0H)2 abundance 
contours show inflections along the CH4-CO boundary, because Mg(0H)2 is sensitive to 
the water abundance in reaction (1221) . Within either the CH4 or CO field, the Mg(0H)2 
abundance is pressure-independent. In principle, the abundance of Mg(0H)2 is therefore 
diagnostic of atmospheric temperature if the object's metallicity is known. 

At low temperatures and high pressures, Mg(0H)2 becomes the most abundant Mg- 
bearing gas, as illustrated in Figure |5l Magnesium hydroxide is thus e xpected to be the 



domi nant magnesium gas below ~1550 K in the atmosphere of Jupiter (cf. iFegley fc Lodders 



19941 ) and below ~980 K on Gliese 229B. 
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4-2.5. Magnesium Monoxide, MgO 

The abundance of MgO as a function of pressure and temperature is illustrated in Figure 
[7^. Magnesium monoxide forms via the net thermochemical reaction 

Mg + H2O = MgO + H2 (29) 

Below the Mg-silicate cloud deck, at pressures less than ~ 10 bar, the MgO abundance 
is effectively pressure-independent and is therefore potentially diagnostic of atmospheric 
temperature if the metallicity is known. 

For similar reasons as for MgOH (see §4.2.3p . the MgO abundance above the magnesium- 
silicate clouds is independent of metallicity and is unaffected by the CH4-CO boundary. For 
example, the Mg abundance decreases with metallicity, as shown by equations and (12^ . 
whereas the H2O abundance increases with metallicity. These effects cancel out in reaction 
fl29|) . with the result that MgO is unaffected by changes in metallicity or the prevailing 
carbon chemistry. 



4.2.6. Magnesium Sulfide, MgS 

Figure Wp displays mole fraction contours for MgS in protosolar composition gas. The 
equilibrium abundance of MgS is governed by the reaction 

Mg + H2S = MgS + H2 (30) 

Near the silicate cloud base, the MgS abundance is ~ 10 ppb in a solar-met allicity gas and 
decreases at higher temperatures and lower pressures as H2S is replaced by monatomic S. 
Above the Mg-silicate clouds, Xugs is proportional to P^^ and is metallicity-independent 
(see Table 12]). Magnesium sulfide is typically among the more abundant Mg-bearing gases 
in brown dwarfs and extrasolar giant planets, and becomes the second most abundant mag- 
nesium gas at temperatures below ~1600 K in the atmospheres of L dwarfs and Pegasi 
planets. 



4.3. Magnesium Chemistry in Substellar Objects 

Magnesium gas chemistry along the pressure-temperature profiles of Jupiter, Gliese 
229B, HD209458b, and an L dwarf (Teg = 1800 K) are illustrated in FiguresEli-IHli. Monatomic 
Mg is the dominant Mg-bearing gas throughout the atmospheres of Pegasi planets and L 
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dwarfs, and in the deep atmospheres of giant planets and T dwarfs. Furthermore, the Mg 
abundance below the Mg-silicate cloud deck is essentially constant and representative of the 
total Mg abundance in the atmosphere. Upon condensation, the abundances of magnesium- 
bearing gases rapidly decrease with decreasing temperature above the magnesium-silicate 
clouds. In objects with lower effective temperatures, a number of other Mg-bearing gases 
become relatively abundant and Mg(0H)2 replaces Mg as the most abundant magnesium gas 
at the ~ 1550 K level on Jupiter and the ~ 980 K level on Gliese 229B. Magnesium hydride 
(MgH) is the second most abundant magnesium gas in the deep atmospheres of substellar 
objects until it is replaced at lower temperatures by Mg(0H)2 and MgOH (in giant planets 
and T dwarfs) or MgS (in L dwarfs and Pegasi planets). Magnesium hydroxide (Mg(0H)2) 
and MgO are potential atmospheric temperature probes, respectively, above and below the 
magnesium-silicate clouds. 



5. Silicon Chemistry 



5.1. Overview of Silicon Chemistry 

An overview of silicon chemistry as a function of pressure and temperature is illustrated 
in Figure |9l The dotted lines labeled Mg2Si04(s,l) and MgSi03(s,l) show the condensation 
temperatures of forsteri te and enstatit e , whi ch together remove nearly all silicon from the 
overlying atmosphere. ICushing et al.l ( 120061 ) found that silicate absorption features near 
10 yum in mid-type L dwarf spectra are consistent with the presence of these magnesium 
silicates, but noted th e possibility of addi tional absorption by quartz (Si02) grains based 



upon the predictions of lHelling et al.l ( 120061 ). In contrast to the models of iHelling et al.l ( 12006 



20081 ). we find that quartz will not condense in the atmospheres of substellar objects unless 
enstatite condensation is suppressed. This is demonstrated in Figure 10, which shows the 
elemental distribution of Si in condensed phases at 1 bar total pressure in a solar-metallicity 
gas with (Figure [TOb ) and without (Figure ITOb) enstatite condensation. As shown in Figure 
[TOl Mg2Si04 formation consumes nearly half of the total Si abundance because the solar 
elemental abundances of Mg and Si are approximately equal. Enstatite formation plausibly 
proceeds via reactions between SiO gas and pre-existing forsterite grains and continues until 
nearly all silicon is consumed. Thus, quartz (Tcond ~ 1550 K) can only form in the absence 
of enstatite (Tcond ~ 1600 K) , because MgSiOs otherwise efficiently removes silicon from the 
gas phase. Even in the absence of gas-grain reactions between SiO and Mg2Si04, the vapor 
pressures of Mg and SiO above forsterite (reaction [T71) remain high enough to drive enstatite 
cloud formation via the net thermochemical reaction (I19p . so that MgSiOs condenses instead 
of Si02. We therefore conclude that Si02 will not condense within the silicate cloud. 
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The most abundant Si-bearing gas over a wide range of pressures and temperatures is 
SiO, until it is replaced at higher pressures by silane, SiH4 (see Figure [9]). The equilibrium 
conversion between SiO and SiH4 is 

SiH4 + H2O = SiO + 3H2. (31) 

The position of the line where SiO and SiH4 have equal abundances A(SiO) = A(SiH4) 
~ |SSigas (where SSigas is the total amount of silicon in the gas) is given by 

log Pt = 4.68 - 3086/T + 0.5[Fe/H]. (32) 

As temperatures decrease, SiH4 is replaced by SiHaF via the net thermochemical reaction 

SiH4 + HF = SiHsF + H2. (33) 

This reaction is independent of pressure and the position of the SiH4-SiH3F boundary is 
given by 

T = 6948/(7.24 - [Fe/H]), (34) 

and occurs at T ~ 960 K in a solar-metallicity gas. Meanwhile, The equilibrium between 
SiO and SiHaF is 

SiO + HF + 2H2 = SiHaF + H2O, (35) 
where the position of the SiO-SiHsF line is given by 

log Pt = 8.30 - 6560/T, (36) 

independent of the metallicity. The SiH4-SiO, SiH4-SiH3F, and SiO-SiHsF equal abundance 
lines intersect to form a "triple point" at T ~ 960 K and Pt ~ 10^'^'' bar in solar-metallicity 
gas, where all three gases have equal abundances [A(SiO) = A(SiH4) = A(SiH3F) ^ |SSigas]. 

At lower temperatures, SiH2F2 replaces SiHsF via the net thermochemical reaction 

SiHsF + HF = SiH2F2 + H2. (37) 

This reaction is also independent of pressure and the SiH3F-SiH2F2 boundary is given by 

T = 7040/(7.68 - [Fe/H]), (38) 

and is located at T ~ 917 K in a solar-metallicity gas. The conversion between SiO and 
SiH2F2 takes place by the reaction 

SiO + HF + H2 = SiH2F2 + H2O, (39) 
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where the position of the SiO-SiH2F2 boundary is given by 

log Pt = 12.14 - 10080/T - 0.5[Fe/H]. (40) 

In a solar-metalhcity gas, equations fl36|) and f l40|) intersect to form the SiO-SiH3F-SiH2F2 
"triple point" at T ~ 917 K and Pt ~ 10^'^^ bar, where all three gases have equal abundances 
{X ^ 10~^^'^ for [Fe/H] = 0). However, the abundances of SiHaF and SiH2F2 in this region 
are extremely low because most silicon is removed from the atmosphere by cloud formation 
at deeper levels. 



5.2. Chemical Behavior of Silicon-bearing Gases 

5.2.1. Silicon Monoxide, SiO 

The chemical behavior of SiO in a protosolar composition gas is illustrated in Figure [11] 
as a function of pressure and temperature. Within the SiO field, the SiO abundance below 
the magnesium-silicate clouds is given by 

logXsio^ -4.20+ [Fe/H], (41) 

and SiO contains ~100% of the atmospheric silicon inventory (Xsio ~ ^ssi)- Upon Mg- 
silicate condensation, silicon is efficiently removed from the gas phase and the SiO abundance 
rapidly decreases with decreasing temperature. Above the clouds, the SiO abundance is 
governed by its vapor pressure over rock. Curvature in the SiO contour lines along the 
CH4-CO equal abundance boundary results from the effect of carbon chemistry on the H2O 
abundance in reactions ( !T7|) and ( !T9|) . The SiO abundance is thus 

Xsio ^ 13.21 - 28817/T - log Pt - [Fe/H], (42) 

in CH4-dominated objects and 

Xsio ^ 13.67 - 28817/T - log Pt - [Fe/H], (43) 

in CO-dominated objects, inversely proportional to pressure and metallicity. These expres- 
sions, along with chemical equilibria, are used to determine the equilibrium abundances of 
other Si-bearing gases. Expressions giving the abundances of silicon species as a function 
of temperature, pressure, and metallicity in a protosolar composition gas are listed in Table 
El We note that the types of cloud condensates present will affect the gas chemistry of 
SiO and subsequently all other Si-bearing gases. For example, as shown in Figure [121 SiO 
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mole fraction abundances above the clouds are ~0.5 dex higher if enstatite formation is sup- 
pressed and replaced by Si02 condensation (see §5.1l above). The SiO abundance is therefore 
potentially diagnostic of weather and cloud composition in brown dwarf atmospheres. 

Silicon monoxide has not yet been detected in the atmospheres of brown dwarfs or giant 
planets, but has been observed in numerous objects including mol ecular clouds, circumst ellar 
envelopes, the photospheres of late-type stars, and sunspots (e.g.. ICampbell et al.lll995l and 
references therein). Abu ndant SiO gas wa s recently detected in the circumstellar disk of the 
13 Pic analog HD172555 Jlisse et al.ll2"o09[ ). 



5.2.2. Silicon Monosulfide, SiS 

Mole fraction contours for SiS are shown in Figure [T3k for a solar- metallicity gas. Silicon 
monosulfide is formed by the reaction between SiO and H2S: 

SiO + H2S = SiS + H2O (44) 

Below the Mg-silicate cloud deck, Xsis ~ 5-10 ppm in a solar-met allicity gas. The SiS 
abundance decreases a t higher temperature s and lower pressures as H2S is replaced by SH 



and monatomic S (e.g. IVisscher et al.l 120061 ). 



Above the clouds, the SiS abundance is inversely proportional to metallicity and total 
pressure (see Table[3]). The SiS abundance contours in Figure [T^ also display curvature along 
the CH4-CO boundary. This shift is more pronounced for SiS than for SiO because the SiS 
abundance in reaction (14^ depends on the SiO and H2O abundances, both of which are af- 
fected by carbon chemistry. For example, when reaction f H^ is at equilibrium, LeChatelier's 
principle shows that more SiS is produced either by adding SiO (or H2S) or removing H2O. 
When moving from the CH4 field to the CO field in a protosolar composition gas, the SiO 
abundance increases and the H2O abundance decreases, yielding a correspondingly large 
increase in the SiS abundance. Silicon sulfide is expected to be the second most abundant 
silicon-bearing gas (after SiO) in the atmospheres of brown dwarfs and extrasolar giant plan- 
ets. Furthermore, it is a potential tracer of weather in these objects because its maximum 
abundance is typically achieved near the magnesium-silicate cloud base. 
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5.2.3. Monatomic Silicon, Si 

Figure [T3b illustrates the chemical behavior of monatomic Si gas as a function of pressure 
and temperature. The abundance of Si is governed by the net thermochemical reaction 

SiO + H2 = Si + H2O. (45) 

Below the Mg-silicate cloud deck and at pressures less than ~ 10 bar, the Si abundance is 
effectively independent of pressure and metallicity and thus potentially diagnostic of atmo- 
spheric temperature. Above the clouds, the Si abundance rapidly decreases with decreasing 
temperature and is proportional to Pj^^ and (see Table [3]). A number of Si absorption 
bands are observed in the near infrared spectra of lo w-mass dwarf stars . These features 



generally weaken and disappear in late- type M dwarfs ( ICushing et al.ll2005l ). consistent with 



the removal of silicon into Mg2Si04 and MgSiOs clouds deeper in the atmosphere. 



5.2.4. Silylidyne, SiH 

The abundance of SiH as a function of pressure and temperature is illustrated in Figure 
IT^ . The SiH abundance is governed by the net thermochemical reaction 

SiO + I.5H2 = SiH + H2O (46) 

Below the Mg-silicate clouds, the SiH abundance is proportional to P^'^ and is independent 
of metallicity. Above the clouds, the SiH abundance is proportional to P^^'^ and 
(see Table [3]). When moving from the CH4 to the CO field, the SiO abundance increases 
and the H2O abundance decreases in reaction fH^ . Both effects serve to increase the SiH 
abundance, resulting in a shift in the SiH contour lines along the CH4-CO boundary. The 
dependence on metallicity means that for every [Fe/H] = +1 dex increase in metallicity, 
the SiH abundance decreases by a factor of 100 at a given pressure and temperature. 



5.2.5. Silylene, SiH2 

Mole fraction contours for SiH2 are shown in Figure [TDd. The equilibrium between SiO 
and SiH2 is given by the reaction 

SiO + 2H2 = SiHa + H2O (47) 

Beneath the magnesium-silicate cloud deck, the metallicity-independent silylene abundance 
increases with total pressure. Upon rock condensation, the amount of SiH2 remaining in 
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the gas rapidly decreases with decreasing temperature. As for SiH, the SiH2 abundance 
is sensitive to the SiO and H2O abundances in reaction and thus a shift occurs in the 
SiH2 contour hues when moving between the CH4 and CO fields. Within each field, the 
SiH2 abundance is pressure-independent (see Figure [H] and Table [3]), and thus potentially 
diagnostic of temperature if the metallicity is known and SiH2 is thermochemical in origin. 

5.2.6. Silyl Radical, SiH^ 

The chemical behavior of the SiHs radical in a solar-metallicity gas is illustrated in 
Figure [T3t. The abundance of SiHs is governed by the reaction 

SiO + 2.5H2 = SiHg + H2O (48) 

Between the H2-H boundary and the Mg-silicate cloud base, the SiHs is proportional to P^-^. 
Above the clouds, the SiHs abundance is proportional to P^^^ and (see Table [3]). There 
is curvature in the SiHs contour lines along the CH4-CO boundary because the prevailing 
carbon chemistry affects both the SiO and H2O abundances in reaction (HHj) 

5.2.7. Silane, SiH^ 

Mole fraction abundance contours for SiH4 are illustrated in Figure [HJi. The equilibrium 
conversion between SiO and SiH4 is represented by reaction (13T|) . and expressions giving the 
SiH4 abundance as a function of pressure, temperature, and metallicity are listed in Table 
[31 Below the clouds, the silane abundance has a strong (P^) dependence on total pressure. 
Upon Mg-silicate cloud formation, there is an inflection in the SiH4 contour lines, and the 
SiH4 abundance above the clouds is proportional to Pt and m~^. Curvature in the SiH4 
contour lines along the CH4-CO boundary results from the effect of carbon chemistry on the 
SiO and H2O abundances in reaction fl3T]) . 

At high pressures and temperatures (> 960 K for [Fe/H] = 0), SiH4 becomes the 
most abundant sihcon-bearing gas. We thus expect SiH4 to be the dominant silicon gas 
below the 1031 K level in the deep atmosphere of Jupiter ([Fe/H] +0.5). However, we 
again emphasize the efficiency with which rock condensation removes silicon from the gas, 
and we expect an abundance of XsiH4 ~ 10^^^'^ at the 1031 K level on Jupiter. At lower 
temperatures, SiH4 is surpassed by SiHsF and SiH2F2. Silane is expected to be the third 
most abundant Si-bearing gas in the deep atmosphere of Gliese 229B. 

We also point out the trend in pressure dependence for the silicon hydrides below 
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(XsiH oc XsiHa oc P^, XsiHg C)C P-^'^, XsiH4 oc P|) and above the silicate clouds 

(XgiH oc Py°'^, XsiH2 oc P^, XsiHg oc P^'^, -^siH4 OC P^), wliicli is evident in the shapes of 
the contour lines in Figure HH The abundance of each of the silicon hydrides is metallicity- 
independent below the Mg-silicate clouds and proportional to above the clouds. 

5.3. Silicon Chemistry in Substellar Atmospheres 

Figures [T5k-[T5li summarize the chemical behavior of silicon gases along the pressure 
temperature profiles of Jupiter, the T dwarf Gliese 229B, the Pegasi planet HD209458b, 
and an L dwarf (Teg = 1800 K). There is a clear trend in silicon chemistry as a function 
of effective temperature. On Jupiter, SiH4 is the dominant Si-bearing gas throughout the 
deep atmosphere, and is replaced by SiHsF and SiH2F2 at lower temperatures (see Figure 
H]). The second most abundant Si-bearing gas is SiO, followed by SiS and a number of other 
silicon gases. In the warmer atmosphere of Gliese 229B, SiO is the dominant Si-bearing 
gas, followed by SiS and SiH4. The relative importance of SiH4 decreases with increasing 
effective temperature, and SiO and SiS are the most important silicon gases throughout 
the atmospheres of Pegasi planets and L dwarfs. The abundances of all the silicon gases 
shown in Figure [15] rapidly decrease with decreasing temperature above the silicate clouds, 
which explains the non-detection of SiH4 and other silicon species in the atmospheres of 
Jupiter and Saturn. Below the cloud base, the important silicon gases SiO and SiS reach 
their maximum abundance and are therefore potential tracers of weather in brown dwarfs 
and Pegasi planets. Silylene (SiH2) and monatomic silicon (Si) are potentially diagnostic of 
atmospheric temperature, respectively, above and below the magnesium silicate clouds. The 
abundant silicon gases SiO and SiS typically achieve their maximum abundance near the 
cloud base and are potential tracers of weather in brown dwarfs and Pegasi planets. 

6. Summary 

The chemical behavior of iron species in substellar atmospheres is strongly affected by 
Fe metal condensation, which efficiently removes most iron from the gas phase. Similarly, 
most magnesium and silicon is removed from the gas by forsterite (Mg2Si04) and enstatite 
(MgSiOs) cloud formation. The equilibrium abundances of Fe-, Mg-, and Si-bearing gases 
rapidly decrease with increasing altitude (and decreasing temperature) above the clouds. 

Monatomic iron is the dominant Fe-bearing gas throughout the atmospheres of L dwarfs 
and Pegasi planets. Other less abundant iron gases become increasingly important in objects 
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with lower effective temperatures, and Fe(0H)2 replaces Fe at low temperatures in T dwarfs 
and giant planets. Magnesium gas chemistry is similar to that of iron. Monatomic Mg is the 
most abundant magnesium gas throughout the atmospheres of L dwarfs and Pegasi planets 
and in the deep atmospheres of giant planets and T dwarfs, where Mg is replaced by Mg(0H)2 
at lower temperatures. A number of Mg-bearing gases become relatively abundant with 
decreasing effective temperature. Silicon monoxide (SiO) is the most abundant Si-bearing 
gas, followed by SiS, throughout the atmospheres of L dwarfs and Pegasi planets and in the 
deep atmospheres of T dwarfs. In objects with lower effective temperatures, a number of 
other silicon gases become increasingly important and SiH4 is the dominant silicon gas in 
the deep atmosphere of Jupiter. At high pressures and low temperatures SiH4 and SiO are 
replaced by SiHsF and/or SiH2F2. 

The abundances of several Fe-, Mg-, and Si-bearing gases are pressure- independent and 
thus, in principle, diagnostic of atmospheric temperature. These include Fe(0H)2, Mg(0H)2, 
and Si below the clouds and FcO, MgO, and SiH2 above the clouds. In addition, a number 
of gases (e.g. Fe, FeH, FeOH, FeS, Mg, MgH, MgOH, MgS, SiO, SiS) may serve as indicators 
of weather since they generally reach their maximum abundance just below the iron metal 
or magnesium-silicate cloud decks. This may be particularly useful for late M dwarfs and 
early L dwarfs in which the metal and silicate clouds are located at relatively shallow depths 
below the photosphere. 
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Table 1. Iron Gas Abundances 



gas M log Xm ~ reaction no. 



Below iron cloud 



Fe 


-4.24+ [Fe/H] 




FeH 


-6.13 - 1905/T + 0.5 log Pt + [Fe/H] 


El 


FeOH 


-9.80 + 2468/T + 0.5 log + 2[Fe/H] 


M 


Fe(0H)2 


-16.64 + 12088/T + logPr + 3[Fe/H] 


m 


FeO 


-7.20 -4713/T + 2[Fe/H] 


m 


FeS 


-8.52 + 964/r + 2[Fe/H] 


m 


Above iron 


cloud 




Fe 


7.23 - 20995/T- log Pt 


m 


FeH 


5.34 - 22900/T- 0.5 log Pt 


m 


FeOH 


5.02 - 18527/T - 0.5 log Pr + 2[Fe/H] + logX^^o 


m 


Fe(0H)2 


1.53 - 8907/r + 2 [Fe/H] + 21ogX^^o 


m 


FeO 


7.62 - 25708/T - log Pt + [Fe/H] + logX(j^o 


m 


FeS 


2.96 - 20031/T - log Pt + [Fe/H] 


m 



Note. — -^RiO defined as the H2O mole fraction in a solar-metallicity 
gas, where logX^^Q ^ —3.12 within the CH4 field and ~ —3.58 within the CO 
field. The effect of metallicity on the atmospheric water abundance (Xu^o) is 
included in the abundance equations. Expressions are valid for temperatures 
between 800 and 2500 K and metallicities up to [Fe/H] = +0.5 dex. 
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Table 2. Magnesium Gas Abundances 



gas M log Xm ~ reaction no. 



Below magnesium-silicate clouds 

Mg -4.15 + [Fe/H] 

MgH -5.46 - 4236/r + 0.5 log Pt + [Fe/H] [27] 

MgOH -8.98 + 1672/T + 0.5 log PT + 2[Fe/H] M 

Mg(0H)2 -15.54 + 10267/T + log Pt + 3[Fe/H] [22] 

MgO -6.12 - 7306/r + 2 [Fe/H] M 

MgS -7.94 + 2 [Fe/H] M 

Above magnesium-silicate clouds 

Mg 8.25 - 27250/r - log Pt - [Fe/H] - log X^j^q [H] M 

MgH 6.94 - 31486/r - 0.5 log Pt - [Fe/H] - logX^j^o [27] 

MgOH 6.75 - 25578/T- 0.5 log Pt M 

Mg(0H)2 3.53 - 16983/T+ [Fe/H] + log X{j^o M 

MgO 9.63 - 34556/T - log Pt M 

MgS 4.40 - 27250/T - log Pt - log X^j^q M 



Note. — ^HaO defined as the H2O mole fraction in a solar-metallicity 
gas, where logX^j^o ~ —3.12 within the CH4 field and ~ —3.58 within the CO 
field. The effect of metallicity on the atmospheric water abundance (XH20) is 
included in the abundance equations. Expressions are valid for temperatures 
between 800 and 2500 K and metallicities up to [Fe/H] = -1-0.5 dex. 
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Table 3. Silicon Gas Abundances 



gas M log Xm ~ reaction no. 



Below magnesium-silicate clouds; within SiO field 

SiO -4.20 + [Fe/H] 

SiS -5.59 + 666/r + 2[Fe/H] M 

Si -0.44 - 44738/T M 

SiH -2.64 - 11500/T + 0.5 log Pt M 

SiH2 -5.83 - 6422/T + log Pt M 

SiHg -8.97- 2770 /T + 1.5 log Pt SSI 

SiH4 -13.33 + 6172/T + 2 log Pr ED 

SiHaF -20.57 + 13120/T + 2 log Pr + [Fe/H] |35] 

SiHaFa -28.25 + 20160/T + 2 log Pt + 2[Fe/H] M 

Above magnesium-silicate clouds 

SiO 10.09 - 28817/T - logPr - [Fe/H] - logX^j^o [13 M 

SiS 5.38 - 28151/T - log Pp - [Fe/H] - 2 logX^j^o M 

Si 10.45 - 44659/T - logP^ - 2 [Fe/H] - 21ogJ\:(i^o SSI 

SiH 8.30 -40317/T-0.51ogPr-2[Fe/H] -21ogX^^o M 

SiHs 5.11 - 35239/T - 2 [Fe/H] - 2\ogX{^^Q M 

SiHg 1.97-31587/T + 0.51ogPr-2[Fe/H] -21ogX^^o M 

SiH4 -2.39 - 22645 /T + log Pr - 2 [Fe/H] - 21ogX^^o ED 

SiHsF -9.63 - 15697/T + log Pr - [Fe/H] - 2 log X^j^q 

SiHsFs -17.31 - 8657/T + log Pr - 2 log X{^^q M 



Note. — o is defined as the H2O mole fraction in a solar-metallicity 
gas, where logX^j^Q ^ —3.12 within the CH4 field and ~ —3.58 within the CO 
field. The effect of metallicity on the atmospheric water abundance (Xnao) is 
included in the abundance equations. Expressions are valid for temperatures 
between 800 and 2500 K and metallicities up to [Fe/H] = +0.5 dex. 



-28- 




log (P^ bar) 



Fig. 1. — Overview of iron chemistry as a function of temperature and pressure in a protosolar 
composition gas. The sohd hne indicates where Fe and Fe(0H)2 have equal abundances. The 
dotted hne labeled Fc(s,l) is the condensation curve for iron, with an open circle denoting its 
normal melting point (1809 K). Also shown arc model atmospheric profiles for representative 
substellar objects (dashed lines). See text for details. 



-29- 




Fig. 2. — Mole fraction contours (on a logarithmic scale) for (a) monatomic iron (Fe), (b) 
iron hydride (FeH), (c) iron monohydroxide (FeOH), and (d) iron hydroxide (Fe(0H)2) as a 
function of pressure and temperature in a solar-metallicity gas. The kinks in the abundance 
contours are due to Fe metal or liquid condensation. 
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Fig. 3. — Mole fraction contours (on a logarithmic scale) for (a) iron monoxide (FeO) and 
(b) iron sulfide (FeS) as a function of pressure and temperature in a solar-metallicity gas. 
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Fig. 4. — Iron chemistry along the atmospheric pressure-temperature profile from 1000 to 
2500 K for (a) Jupiter, (b) Ghese 229B, (c) HD209458b, and (d) an L dwarf {T^s = 1800 K). 
Breaks in the Fe curves show where Fe condenses (b,d). Iron is condensed throughout the 
entire temperature range shown here for Jupiter and HD209458b (a,c). 
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Fig. 5. — Overview of magnesium chemistry as a function of temperature and pressure in 
a protosolar composition gas. The solid line indicates where Mg and Mg(0H)2 have equal 
abundances. The dotted lines labeled Mg2Si04(s,l) and MgSi03(s,l) show the condensation 
temperatures of forsterite and enstatite, with circles denoting their normal melting points of 
1803 K (enstatite) and 2163 K (forsterite). The dotted line labeled MgO(s) is the conden- 
sation curve for periclase. See text for details. 
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Fig. 6. — Mole fraction contours (on a logarithmic scale) for (a) monatomic magnesium (Mg), 
(b) magnesium hydride (MgH), (c) magnesium monohydroxide (MgOH), and (d) magnesium 
hydroxide (Mg(0H)2) as a function of pressure and temperature in a solar-metallicity gas. 
The kinks in the abundance contours occur at the condensation temperature of Mg2Si04 
(forsterite) . 
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Fig. 7. — Mole fraction contours (on a logarithmic scale) for (a) magnesium monoxide 
(MgO) and (b) magnesium sulfide (MgS) as a function of pressure and temperature in a 
solar-metallicity gas. 
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Fig. 8. — Magnesium chemistry along the atmospheric pressure-temperature profile from 
1000 to 2500 K for (a) Jupiter, (b) Ghese 229B, (c) HD209458b, and (d) an L dwarf 
(Teflf = 1800 K). Breaks in the Mg curves show where magnesium-sihcates condense (b,c,d). 
Forsterite is condensed throughout the entire temperature range shown here for Jupiter (a). 
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Fig. 9. — Overview of silicon chemistry as a function of temperature and pressure in a 
protosolar composition gas. The sohd hne indicates where dominant Si-bearing gases have 
equal abundances. The dotted hnes labeled Mg2Si04 and MgSiOa show the condensation 
temperatures of forsterite and enstatite, with circles denoting their melting points. Also 
shown are model atmospheric profiles for representative substellar objects (dashed lines). 
See text for details. 
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Fig. 10. — Elemental distribution of silicon into condensed phases at 1 bar in a solar- 
metallicity gas with enstatite (MgSiOa) formation a) included or b) suppressed. Quartz 
(Si02) condensation will only proceed in the absence of enstatite, which otherwise efficiently 
removes sihcon from the gas phase. The silicates mehhte, anorthite, and diopside together 
consume up to ~12% of the total Si abundance. SiO and SiS are the dominant Si-bearing 
gases before removal by condensation. 



-38- 




Fig. 11. — Mole fraction contours (on a logarithmic scale) for silicon monoxide (SiO) as a 
function of pressure and temperature in a solar-metallicity gas. 
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Fig. 12. — Mole fraction abundance of SiO gas at 1 bar total pressure in a solar-met allicity gas 
with enstatite formation included (black line) or suppressed (gray line) above the forsterite 
cloud layer. The formation of Si02 instead of MgSiOa would yield larger mole fraction 
abundances of SiO gas above the clouds. 
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Fig. 13. — Mole fraction contours (on a logarithmic scale) for (a) silicon monosulfide (SiS) 
and and (b) monatomic silicon (Si) as a function of pressure and temperature in a solar- 
metallicity gas. 
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Fig. 14. — Mole fraction contours (on a logarithmic scale) for (a) silylidyne (SiH), (b) 
silylene (SiH2), (c) the silyl radical (SiHs), and (d) silane (SiH4) as a function of pressure 
and temperature in a solar-metallicity gas. The bends in the abundance contours occur 
where magnesium sihcate clouds (MgSiOs, Mg2Si04) condense. 
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Fig. 15. — Silicon chemistry along the atmospheric pressure-temperature profile from 1000 
to 2500 K for (a) Jupiter, (b) Gliese 229B, (c) HD209458b, and (d) an L dwarf {T^s = 1800 
K). Breaks in the SiO curves show where magnesium-sihcates condense (b,c,d). Forsterite is 
condensed throughout the entire temperature shown here for Jupiter (a). 



